High frequency of transition to transversion ratio in the stem region of RNA secondary structure of untranslated region of SARS-CoV-2.

Introduction: The propensity of nucleotide bases to form pairs, causes folding and the formation of secondary structure in the RNA. Therefore, purine (R): pyrimidine (Y) base-pairing is vital to maintain uniform lateral dimension in RNA secondary structure. Transversions or base substitutions between R and Y bases, are more detrimental to the stability of RNA secondary structure, than transitions derived from substitutions between A and G or C and T. The study of transversion and transition base substitutions is important to understand evolutionary mechanisms of RNA secondary structure in the 5'  and 3'  untranslated (UTR) regions of SARS-CoV-2. In this work, we carried out comparative analysis of transition and transversion base substitutions in the stem and loop regions of RNA secondary structure of SARS-CoV-2. Methods: We have considered the experimentally determined and well documented stem and loop regions of 5' and 3' UTR regions of SARS-CoV-2 for base substitution analysis. The secondary structure comprising of stem and loop regions were visualized using the RNAfold web server. The GISAID repository was used to extract base sequence alignment of the UTR regions. Python scripts were developed for comparative analysis of transversion and transition frequencies in the stem and the loop regions. Results: The results of base substitution analysis revealed a higher transition (ti) to transversion (tv) ratio (ti/tv) in the stem region of UTR of RNA secondary structure of SARS-CoV-2 reported during the early stage of the pandemic. The higher ti/tv ratio in the stem region suggested the influence of secondary structure in selecting the pattern of base substitutions. This differential pattern of ti/tv values between stem and loop regions was not observed among the Delta and Omicron variants that dominated the later stage of the pandemic. It is noteworthy that the ti/tv values in the stem and loop regions were similar among the later dominant Delta and Omicron variant strains which is to be investigated to understand the rapid evolution and global adaptation of SARS-CoV-2. Conclusion: Our findings implicate the lower frequency of transversions than the transitions in the stem regions of UTRs of SARS-CoV-2. The RNA secondary structures are associated with replication, translation, and packaging, further investigations are needed to understand these base substitutions across different variants of SARS-CoV-2.

Machine learning approach identifies prominent codons from different degenerate groups influencing gene expression in bacteria.

Unequal usage of synonymous codons is known as codon usage bias (CUB), which is generally different between the high-expression genes (HEG) and low-expression genes (LEG) in organisms is not yet adequately reported across different bacteria. In this study, a machine learning-based approach was implemented initially to find out codons that are significantly different between the HEG and LEG in Escherichia coli. It identified Cys codons such as UGU and UGC, Lys codons such as AAA and AAG that were least influenced by gene expression. Codons such as UCU (Ser), CUG (Leu), GGG (Gly), CGG (Arg) etc. were identified to be influenced maximum by the gene expression. The study was extended to analyze codon usage in 683 other bacterial species. Cys (UGU/UGC) and Ser (AGU/AGC) codons were identified being the least different between the two groups of genes across these bacterial species. Codons such as CGA, CUG, GGG, GCC, ACC, AUA, and AUC were identified to be influenced by the gene expression across majority of these species. This study supports the role of CUB on gene expression across bacteria and demonstrates a commonality among bacteria regarding behavior of certain codons with regard to gene expression.

Incorporation of transition to transversion ratio and nonsense mutations, improves the estimation of the number of synonymous and non-synonymous sites in codons.

A common approach to estimate the strength and direction of selection acting on protein coding sequences is to calculate the dN/dS ratio. The method to calculate dN/dS has been widely used by many researchers and many critical reviews have been made on its application after the proposition by Nei and Gojobori in 1986. However, the method is still evolving considering the non-uniform substitution rates and pretermination codons. In our study of SNPs in 586 genes across 156 Escherichia coli strains, synonymous polymorphism in 2-fold degenerate codons were higher in comparison to that in 4-fold degenerate codons, which could be attributed to the difference between transition (Ti) and transversion (Tv) substitution rates where the average rate of a transition is four times more than that of a transversion in general. We considered both the Ti/Tv ratio, and nonsense mutation in pretermination codons, to improve estimates of synonymous (S) and non-synonymous (NS) sites. The accuracy of estimating dN/dS has been improved by considering the Ti/Tv ratio and nonsense substitutions in pretermination codons. We showed that applying the modified approach based on Ti/Tv ratio and pretermination codons results in higher values of dN/dS in 29 common genes of equal reading-frames between E. coli and Salmonella enterica. This study emphasizes the robustness of amino acid composition with varying codon degeneracy, as well as the pretermination codons when calculating dN/dS values.

Stem Region of tRNA Genes Favors Transition Substitution Towards Keto Bases in Bacteria.

Transversion and transition mutations have variable effects on the stability of RNA secondary structure considering that the former destabilizes the double helix geometry to a greater extent by introducing purine:purine (R:R) or pyrimidine:pyrimidine (Y:Y) base pairs. Therefore, transversion frequency is likely to be lower than that of transition in the secondary structure regions of RNA genes. Here, we performed an analysis of transition and transversion frequencies in tRNA genes defined well with secondary structure and compared with the intergenic regions in five bacterial species namely Escherichia coli, Klebsiella pneumoniae, Salmonella enterica, Staphylococcus aureus and Streptococcus pneumoniae using a large genome sequence data set. In general, the transversion frequency was observed to be lower than that of transition in both tRNA genes and intergenic regions. The transition to transversion ratio was observed to be greater in tRNA genes than that in the intergenic regions in all the five bacteria that we studied. Interestingly, the intraspecies base substitution analysis in tRNA genes revealed that non-compensatory substitutions were more frequent than compensatory substitutions in the stem region. Further, transition to transversion ratio in the loop region was observed to be significantly lesser than that among the non-compensatory substitutions in the stem region. This indicated that the transversion is more deleterious than transition in the stem regions. In addition, substitutions from amino bases (A/C) to keto bases (G/T) were also observed to be more than the reverse substitutions in the stem region. Substitution from amino bases to keto bases are likely to facilitate the stable G:U pairing unlike the reverse substitution that facilitates the unstable A:C pairing in the stem region of tRNA. This work provides additional support that the secondary structure of tRNA molecule is what drives the different substitutions in its gene sequence.

The Entner-Doudoroff and Nonoxidative Pentose Phosphate Pathways Bypass Glycolysis and the Oxidative Pentose Phosphate Pathway in Ralstonia solanacearum.

In Ralstonia solanacearum, a devastating phytopathogen whose metabolism is poorly understood, we observed that the Entner-Doudoroff (ED) pathway and nonoxidative pentose phosphate pathway (non-OxPPP) bypass glycolysis and OxPPP under glucose oxidation. Evidence derived from 13C stable isotope feeding and genome annotation-based comparative metabolic network analysis supported the observations. Comparative metabolic network analysis derived from the currently available 53 annotated R. solanacearum strains, including a recently reported strain (F1C1), representing the four phylotypes, confirmed the lack of key genes coding for phosphofructokinase (pfk-1) and phosphogluconate dehydrogenase (gnd) enzymes that are relevant for glycolysis and OxPPP, respectively. R. solanacearum F1C1 cells fed with [13C]glucose (99% [1-13C]glucose or 99% [1,2-13C]glucose or 40% [13C6]glucose) followed by gas chromatography-mass spectrometry (GC-MS)-based labeling analysis of fragments from amino acids, glycerol, and ribose provided clear evidence that rather than glycolysis and the OxPPP, the ED pathway and non-OxPPP are the main routes sustaining metabolism in R. solanacearum The 13C incorporation in the mass ions of alanine (m/z 260 and m/z 232), valine (m/z 288 and m/z 260), glycine (m/z 218), serine (m/z 390 and m/z 362), histidine (m/z 440 and m/z 412), tyrosine (m/z 466 and m/z 438), phenylalanine (m/z 336 and m/z 308), glycerol (m/z 377), and ribose (m/z 160) mapped the pathways supporting the observations. The outcomes help better define the central carbon metabolic network of R. solanacearum that can be integrated with 13C metabolic flux analysis as well as flux balance analysis studies for defining the metabolic phenotypes.IMPORTANCE Understanding the metabolic versatility of Ralstonia solanacearum is important, as it regulates the trade-off between virulence and metabolism (1, 2) in a wide range of plant hosts. Due to a lack of clear evidence until this work, several published research papers reported on the potential roles of glycolysis and the oxidative pentose phosphate pathway (OxPPP) in R. solanacearum (3, 4). This work provided evidence from 13C stable isotope feeding and genome annotation-based comparative metabolic network analysis that the Entner-Doudoroff pathway and non-OxPPP bypass glycolysis and OxPPP during the oxidation of glucose, a component of the host xylem pool that serves as a potential carbon source (5). The outcomes help better define the central carbon metabolic network of R. solanacearum that can be integrated with 13C metabolic flux analysis as well as flux balance analysis studies for defining the metabolic phenotypes. The study highlights the need to critically examine phytopathogens whose metabolism is poorly understood.

Codon degeneracy and amino acid abundance influence the measures of codon usage bias: improved Nc (N̂c ) and ENCprime (N̂'c ) measures.

Effective number of codons (N^c) and its variant N^'c (effective number of codons prime) are the two widely used methods for measuring unequal usage of synonymous codons in coding sequences, known as the codon usage bias (CUB). The mathematical formula used in calculating N^c and N^'c values is giving inappropriate measures of CUB in case of low abundance of amino acids. In addition, the magnitude of error also varies according to codon degeneracy. In this study, a modified formula for N^c and N^'c has been developed to measure the CUB more accurately. Online implementations of the modified formula are available in the web portal at http://agnigarh.tezu.ernet.in/~ssankar/cub.php.

Comparative analysis of codon usage bias in Crenarchaea and Euryarchaea genome reveals differential preference of synonymous codons to encode highly expressed ribosomal and RNA polymerase proteins.

The present study was undertaken to investigate the pattern of optimal codon usage in Archaea. Comparative analysis was executed to understand the pattern of codon usage bias between the high expression genes (HEG) and the whole genomes in two Archaeal phyla, Crenarchaea and Euryarchaea. The G+C% of the HEG was found to be less in comparison to the genome G+C% in Crenarchaea, whereas reverse was the case in Euryarchaea. The preponderance of U/A ending codons that code for HEG in Crenarchaea was in sharp contrast to the C/G ended ones in Euryarchaea. The analysis revealed prevalence of Uending codons even within theWWY(nucleotide ambiguity code) families in Crenarchaea vis-à-vis Euryarchaea, bacteria and Eukarya. No plausible interpretation of the observed disparity could be made either in the context of tRNA gene composition or genome G+C%. The results in this study attested that the preferential biasness for codons in HEG of Crenarchaea might be different from Euryarchaea. The main highlights are (i) varied CUB in the HEG and in the whole genomes in Euryarchaea and Crenarchaea. (ii) Crenarchaea was found to have some unusual optimal codons (OCs) compared to other organisms. (iii) G+C% (and GC3) of the HEG were different from the genome G+C% in the two phyla. (iv) Genome G+C% and tRNA gene number failed to explain CUB in Crenarchaea. (v) Translational selection is possibly responsible for A+T rich OCs in Crenarchaea.

Discrepancy among the synonymous codons with respect to their selection as optimal codon in bacteria.

The different triplets encoding the same amino acid, termed as synonymous codons, are not equally abundant in a genome. Factors such as G + C% and tRNA are known to influence their abundance in a genome. However, the order of the nucleotide in each codon per se might also be another factor impacting on its abundance values. Of the synonymous codons for specific amino acids, some are preferentially used in the high expression genes that are referred to as the 'optimal codons' (OCs). In this study, we compared OCs of the 18 amino acids in 221 species of bacteria. It is observed that there is amino acid specific influence for the selection of OCs. There is also influence of phylogeny in the choice of OCs for some amino acids such as Glu, Gln, Lys and Leu. The phenomenon of codon bias is also supported by the comparative studies of the abundance values of the synonymous codons with same G + C. It is likely that the order of the nucleotides in the triplet codon is also perhaps involved in the phenomenon of codon usage bias in organisms.

Selection on GGU and CGU codons in the high expression genes in bacteria.

The fourfold degenerate site (FDS) in coding sequences is important for studying the effect of any selection pressure on codon usage bias (CUB) because nucleotide substitution per se is not under any such pressure at the site due to the unaltered amino acid sequence in a protein. We estimated the frequency variation of nucleotides at the FDS across the eight family boxes (FBs) defined as Um(g), the unevenness measure of a gene g. The study was made in 545 species of bacteria. In many bacteria, the Um(g) correlated strongly with Nc'-a measure of the CUB. Analysis of the strongly correlated bacteria revealed that the U-ending codons (GGU, CGU) were preferred to the G-ending codons (GGG, CGG) in Gly and Arg FBs even in the genomes with G+C % higher than 65.0. Further evidence suggested that these codons can be used as a good indicator of selection pressure on CUB in genomes with higher G+C %.

Constraint on di-nucleotides by codon usage bias in bacterial genomes.

It has been reported earlier that the relative di-nucleotide frequency (RDF) in different parts of a genome is similar while the frequency is variable among different genomes. So RDF is termed as genome signature in bacteria. It is not known if the constancy in RDF is governed by genome wide mutational bias or by selection. Here we did comparative analysis of RDF between the inter-genic and the coding sequences in seventeen bacterial genomes, whose gene expression data was available. The constraint on di-nucleotides was found to be higher in the coding sequences than that in the inter-genic regions and the constraint at the 2nd codon position was more than that in the 3rd position within a genome. Further analysis revealed that the constraint on di-nucleotides at the 2nd codon position is greater in the high expression genes (HEG) than that in the whole genomes as well as in the low expression genes (LEG). We analyzed RDF at the 2nd and the 3rd codon positions in simulated coding sequences that were computationally generated by keeping the codon usage bias (CUB) according to genome G+C composition and the sequence of amino acids unaltered. In the simulated coding sequences, the constraint observed was significantly low and no significant difference was observed between the HEG and the LEG in terms of di-nucleotide constraint. This indicated that the greater constraint on di-nucleotides in the HEG was due to the stronger selection on CUB in these genes in comparison to the LEG within a genome. Further, we did comparative analyses of the RDF in the HEG rpoB and rpoC of 199 bacteria, which revealed a common pattern of constraints on di-nucleotides at the 2nd codon position across these bacteria. To validate the role of CUB on di-nucleotide constraint, we analyzed RDF at the 2nd and the 3rd codon positions in simulated rpoB/rpoC sequences. The analysis revealed that selection on CUB is an important attribute for the constraint on di-nucleotides at these positions in bacterial genomes. We believe that this study has come with major findings of the role of CUB on di-nucleotide constraint in bacterial genomes.

Transfer RNA gene numbers may not be completely responsible for the codon usage bias in asparagine, isoleucine, phenylalanine, and tyrosine in the high expression genes in bacteria.

It is generally believed that the effect of translational selection on codon usage bias is related to the number of transfer RNA genes in bacteria, which is more with respect to the high expression genes than the whole genome. Keeping this in the background, we analyzed codon usage bias with respect to asparagine, isoleucine, phenylalanine, and tyrosine amino acids. Analysis was done in seventeen bacteria with the available gene expression data and information about the tRNA gene number. In most of the bacteria, it was observed that codon usage bias and tRNA gene number were not in agreement, which was unexpected. We extended the study further to 199 bacteria, limiting to the codon usage bias in the two highly expressed genes rpoB and rpoC which encode the RNA polymerase subunits ß and ß', respectively. In concordance with the result in the high expression genes, codon usage bias in rpoB and rpoC genes was also found to not be in agreement with tRNA gene number in many of these bacteria. Our study indicates that tRNA gene numbers may not be the sole determining factor for translational selection of codon usage bias in bacterial genomes.

Higher tRNA diversity in thermophilic bacteria: a possible adaptation to growth at high temperature.

We have done a comparative study of tRNA diversity and total tRNA genes among different strains of bacteria with respect to the optimum growth temperature of the cells. Our observation suggests that higher tRNA diversity usually occurs in thermophiles in comparison to non-thermophiles. Among psychrophiles total tRNA was observed to be more than two-fold higher than in the non-psychrophiles. Though tRNA diversity and total tRNA have recently been shown to be affected by an organism's genomic GC% and growth rate, this work is the first report on growth temperature affecting these features in bacteria. This work extends the list of molecular features undergoing adaptation due to growth temperature and supports the view that growth temperature acts as a selecting factor at the molecular level during evolution.

Variable correlation of genome GC% with transfer RNA number as well as with transfer RNA diversity among bacterial groups: alpha-proteobacteria and tenericutes exhibit strong positive correlation.

Mainly two selective forces are thought to affect transfer RNA (tRNA) gene numbers in cells: firstly, transfer RNA genes have coevolved with codon usage to optimize the translation rate in organisms; and secondly, tRNA gene number is positively correlated with growth rate of organisms. Since a codon ending with 'G' is recognized with lower efficiency than the synonymous codon ending with 'A' by a tRNA, organisms with high genome GC% are expected to possess tRNA molecules with 'C' at the first anticodon position to decode efficiently the former codon. This is in concordance with the observation of higher tRNA diversity in several bacteria with higher genome GC%. In addition to this, tRNA gene redundancy has been reported to be correlated with genome size. In this paper we have carried out a comprehensive correlation analysis of tRNA gene diversity, total tRNA number, genome GC% and genome size, in 626 bacterial strains (belonging to 433 species). A common positive trend of correlation in bacteria of GC% with tRNA diversity, of GC% with genome size and of genome size with total tRNA number, is observed, even within individual bacterial groups. Comparison of organisms' growth rate data with their genome GC% indicated that organisms with higher growth rate are most likely with <50% GC. This work extends the list of features in organisms that relate to genome GC%.

A study in entire chromosomes of violations of the intra-strand parity of complementary nucleotides (Chargaff's second parity rule).

Chargaff's rule of intra-strand parity (ISP) between complementary mono/oligonucleotides in chromosomes is well established in the scientific literature. Although a large numbers of papers have been published citing works and discussions on ISP in the genomic era, scientists are yet to find all the factors responsible for such a universal phenomenon in the chromosomes. In the present work, we have tried to address the issue from a new perspective, which is a parallel feature to ISP. The compositional abundance values of mono/oligonucleotides were determined in all non-overlapping sub-chromosomal regions of specific size. Also the frequency distributions of the mono/oligonucleotides among the regions were compared using the Kolmogorov-Smirnov test. Interestingly, the frequency distributions between the complementary mono/oligonucleotides revealed statistical similarity, which we named as intra-strand frequency distribution parity (ISFDP). ISFDP was observed as a general feature in chromosomes of bacteria, archaea and eukaryotes. Violation of ISFDP was also observed in several chromosomes. Chromosomes of different strains belonging a species in bacteria/archaea (Haemophilus influenza, Xylella fastidiosa etc.) and chromosomes of a eukaryote are found to be different among each other with respect to ISFDP violation. ISFDP correlates weakly with ISP in chromosomes suggesting that the latter one is not entirely responsible for the former. Asymmetry of replication topography and composition of forward-encoded sequences between the strands in chromosomes are found to be insufficient to explain the ISFDP feature in all chromosomes. This suggests that multiple factors in chromosomes are responsible for establishing ISFDP.